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ABSTRACT 



Context. In the current paradigm of cosmic structure formation, dark matter plays a key role on the formation and evolution of galaxies 
through its gravitational influence. On microscopic scales, dark matter particles are expected to annihilate amongst themselves into 
different products, with some fraction of the energy being transferred to the baryonic component. 

Aims. It is the aim of the present work to show that, in the innermost regions of dark matter halos, heating by dark matter annihilation 
may be comparable to the cooling rate of the gas. 

Methods. We use analytical models of the dark matter and gas distributions in order to estimate the heating and cooling rates, as well 
as the energy available from supernova explosions. 

Results. Depending on the model parameters and the precise nature of dark matter particles, the injected energy may be enough 
to balance radiative cooling in the cores of galaxy clusters. On galactic scales, it would inhibit star formation more efficiently than 
supernova feedback. 

Conclusions. Our results suggest that dark matter annihilation prevents gas cooling and star formation within at least 0.01 - 1 per 
cent of the virial radius. 
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1. Introduction 

One of the most remarkable achievements of modern 
Cosmology is the measurement of t he fundamental cons tituents 
of the Universe. Over 80 per cent dSpergel et alj 120060 of the 
matter (one fifth of the total energy density) is currently be- 
lieve d to be composed of non-baryonic dark matter particles (see 
e.g. iBertone et alj |2005, for a review of candidates). The cold 
dark matter scenario has been extremely successful in explain- 
ing many of the observed properties of galaxies over a broad 
range of scales and environments, but there are nevertheless sev- 
eral issues that still defy our understanding. 

One of them is the number of dwarf galaxies orbiting around 
the Milky Way, with numerical simulations predicti ng one or two 
order s of magnitude mor e satellites than observed dKlvpin et alj 
1999; Mo ore et all 1999). Perhaps the currently mo st favored 
explanation is that photoionization dEfstathioull 19921) prevented 
gas co oling in the smallest objects. According to lKravtsov et alj 
( 2004), star formation should be strongly suppressed for all halos 
smaller than 10 9 M Q . A similar problem might exist on galaxy 
cluster scal es, which would push the threshold to 10 11 M Q or 
even larger dKase et a l. 2006). 

Actually, observations of the conditional luminosity func- 
tion indicate that the mass-to-light ratio reaches a minimum for 
halo masses around 3 x 10 11 M Q , with objects below 10 10 M 
virtually devoid of galaxies dvan den Bosch et ail 12003). It is 
at present unclear whether the cosmic ultraviolet background 
could provide enough photons to achieve such effect; numeri- 
cal experiments suggest that some fraction of the gas is still ex- 
pected to cool and collapse into stars at the centers of most halos 
dHoeft et al.ir2005l) . 

There is an observed upper mass limit above which star for- 
mation seems to be suppressed as well. This threshold increases 



with redshift (see e.g. iBundv et alj 120061) . in apparent contra- 
diction with the hierarchical picture. Star formation in the most 
massive galaxies takes place at relatively early times and then it 
suddenly shuts off, while less massive objects tend to form their 
stars at later times dCowie et alj|1996l) . 

The problem is particularly noticeable in galaxy clusters, 
where it is difficult to understand why the gas is currently cool- 
ing a t a much slower rat e than expected from its X-ray luminos- 
ity dPeterson et alj 2001). and an external heat source, such as an 
active galactic nucleus, is often invoked in order to explain the 
phenomenon. 

A similar mechanism may also be responsible for the red 
and blue sequences observed in the color-magnitude diagram of 
galaxies, that can only be reproduced by qu enching star forma- 
tion in halos more massive than ~ 10 12 M dCroton et al.ll2006h 
ICattaneo et"alll2006l) . 

Here we propose that annihilation of dark matter particles 
may provide a considerable amount of energy, which, if trans- 
ferred to the surrounding gas, could help alleviating some of the 
discrepancies outlined above. 

On cluster scales, it has been argued that neutra lino annihi- 
lation may play a role in the cooling flow problem (Totani 2004; 
IColafrancesco et aT]|2006l) . and there have been several studies 
assessing the impact of dark matt er decay and/or annihilation 
on th e cosmic ionization h i story ( Padmanabha n & Finkbeinerl 
120051; iMapelli et all 120061; IZhang et all 120061) and the soft 
gamm a-ray background dAhn & Komatsul 120051 ; iRasera et alj 
2006 ), as well as o n the first generation of galaxies 
( Ripa monti et~aT]|2006bl) . 

It is our aim to show that, for a relatively broad range of sce- 
narios, dark matter annihilation may also influence galaxy for- 
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mation and evolution by quenching gas cooling and star forma- 
tion near the center of dark matter halos. 



2. Heat from cold dark matter 

When two dark matter particles annihilate, all their energy goes 
into the annihilation products. Some fraction, f m i, will be radi- 
ated away at different wavelengths, from gamma rays to radio, 
and it c an be used to impo se constraints on the nature of dark 
matter dBertone et al.l 12005b . The present study focuses on the 
remaining fraction, / a b s = 1 - / ra d, that is eventually absorbed by 
the surrounding baryonic gas, which is thus heated at a rate 



; /abs (crv) 2m dm c 



/abs Cp 2 dm (crv) c 2 / (2m dm ) ( 1 ) 



where the dot denotes derivative with respect to time, u is the gas 
energy per unit volume, m dm , p dm , and n dm are the mass, density, 
and number density of dark matter particles, respectively, (crv) 
is their annihilation cross-section, and we have assumed in the 
last step non-self-conjugate particles, n dm = "dm* = Pdm/(2»idm)- 

The clumping factor, C = \p 2 , ) / (pdm), accounts for the 
presence of substructure s (Bergst rom et alJI 19991) . Current sim- 
ulations (Die mand et al J 120061) indicate that C should be of 
the order of a factor of 2 or 3, although higher values (e.g. 
IColafrancesco et aTll2006l) . C > 10, cannot be excluded. It is im- 
portant to note, though, that the distribution of subhalos is less 
concentrated than that of the smooth component of the main halo 
( Naga i & Krav"tsovll2005l) . and therefore C should actually be a 
decreasing function of pd m . 

It has been recently shown that, for light dark matter parti- 
cles (m dm ~ 3 — 10 MeV) annihilating in a neutra l unmagnetized 
gas at the average cosmic density, / a b S > 0.03 dRipamonti et alj 
l2006al) . At typical galactic densities, the absorbed fraction is 
expected to increase due to the higher rate of Coulomb colli- 
sions, up to / a bs ~ 1 . However, for heavier candidates (m dm > 
100 MeV) more radiation is expected through inverse Compton 
and synchrotron emission from the initially relativistic annihila- 
tion products. Thermalization becomes again efficient when the 
energy of these particles (most notably electrons and positrons, 
but also protons and other particles) drop s below a few GeV, 
leading to / abs ~ 0.1 for m dm < 100 GeV dTotanill2004l) . 

For thermal relics, the observed cosmic density imposes 
(crv) 3 x 10~ 26 cm 3 s _1 at the time of decoupling. However, 
light particles with 0.511 < m dm < 100 MeV and such a cross- 
section would produce too many positrons and gamma rays 
within our galaxy. Compatibility with IN TEGRAL/SPI mea- 
surements requires (Ascasib ar et alj|2006l) that the present-day 
annihilation cross-section into electron-positron pairs satisfies 



(crvW < 3 x 10 
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Substituting this expression in equation ((T), the energy input 
increases proportionally to m dm until it reaches a maximum at 
fidm ~ 100 MeV, and then it declines as mi. 

For the density profile of the dark matter halo, we adopt the 
general formula 



Pdn 



(r/r s ) a (1 + r/r.) 



3-Q 



(3) 



where p s and r s are the characteristic density and radius of 
the object, and a is the asymptotic logar ithmic slope at the 
center. Cosmological N-body simulations (Nav arro et al.lll99"7h 
suggest a ~ 1. Lower values have been inferred from the 



rotati on curves of dwarf and low surface brightness galax- 
ies dFlores & Primackl 1 1994 iMoord Q994), although it is at 
present unclear whe ther these observations may be consisten t 
with steeper profiles (Havas hi et al.l200 4; Spekk ens et alj2 005). 
On the other ha nd, a > 1 is expected from adiabatic contrac- 
tion due to stars dBlumenthal et alJI 19861) and/or a supermassive 
black hole (Gondo lo & Silkll 199911 Finally, the density profile 



Pdm - Pie 



(4) 



with (3 ^0.18 also d escribes the results of numerical simulations 
dNavarro et al . 2004), and it tends to a finite value at the center. 
In this expression, p2 is the density at the radius r2 where the 
logarithmic slope is equal to -2. For the profile (0, this occurs 
at n = (2 - a)r s . In order to model halos o n different scales, we 
account for the mass-concentration relation (Bullock e t alj200ll) 
according to 
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M200 



20 kpc \ 10 12 M 



0.46 



(5) 



where M200 is the mass enclosed at an overdensity of 200 times 
the critical density. 



3. Results 

We compare in Figure Q] the heating rate given by equation ([]]) 
with the local cooling rate of the gas and the energy available 
from supernova feedback. In order to obtain an upper limit (Case 
I), we set m dm = 100 MeV, (crv) = 3 x 10~ 26 cm 3 s -1 , / abs = 1, 
and C = 10. More realistic values, / a b s =0.1 and C = 1, are 
assumed for Case II (m dm = 10 MeV, (crv) = 3 x 10~ 28 cm 3 s -1 ) 
and Case III (m dm = 100 GeV, (crv) = 3 x 10~ 26 cm 3 s" 1 ). 

The dark matter density profile follows expression ©, and 
the baryonic compo nent is modeled as a p olytrope with effec- 
tive index y =; 1.18 (Ascasibar et al. 2003). The cooling rate of 
the gas is determined by the functi on A(r gas ) = w/p 2 as tabu- 
lated in Sutherland & Dopital (119931) . and new stars are formed 
as dKennicuttll 19981) 



0.02: 



T"dyn 



(6) 



where T dyn = y , (37r)/(16Gp gas ) is the local dynamical timefl 
Supernova explosions can inj ect energy at a r ate u - esNp*, with 
£sn - 4x 10 48 erg M ' for a lSalpeteif dl955h initial mass func- 
tion. 

Heat from annihilating dark matter particles is usually sev- 
eral orders of magnitude below the radiative cooling rate of the 
gas, except in the central part, where the dark matter density 
becomes much larger than the gas density. In this region, dark 
matter annihilation not only provides more energy than super- 
novae; in fact, it would prevent gas cooling and star formation 
completely. 

The heated gas would expand and rise buoyantly, creating 
winds, shocks, and turbulence, and a full three-dimensional sim- 
ulation would be required in order to evaluate the net effect on 
the ambient baryonic medium. Recent numerical studies of the 
effect of cosmic rays accelerated in structure formation shocks 
show that the total mass-to-light ratio of small halos and the faint 



1 If the cooling time r c = u/ii < Td yn , the gas can cool to tempera- 
tures ~ 10 4 K in pressure equilibrium with the ambient medium, and its 
density will be enhanced by a factor r gas /(10 4 K). 
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Fig. 1. Heat injection (dashed lines) for different dark matter 
scenarios (see text), compared to the cooling rate of the gas 
(solid lines) and feedback from supernovae (dotted lines) in ha- 
los of different mass. 



end of the luminosity function can in deed be strongly affe cted by 
the injection of relativistic particles (Jubelg as et alj|2006l) . 

The extent of the heating-dominated region is very sensitive 
to the specific dark matter candidate, which sets the normaliza- 
tion of equation ([]]) through the product of / a bs, (ot), and m^ m . 
We show in Figure [2] the radius r\, at which heating balances 
cooling, as well as the integrated energy injection within r\, and 
rioo, for cases I, II, and III. We also plot the results for a dark 
matter density profile of the form (0 for 0.25 < a < 1 .75. 

When a = 0, dark matter annihilation is not able to coun- 
teract gas cooling at any radius, not even in the most optimistic 
Case I. However, this does not apply to any "cored" density pro- 
file. Our results for expression © are similar to those obtained 
for a = 1 in cases I and II. In Case III, dark matter heating is so 
close to gas cooling that the asymptotic behaviour of the density 
profile has a critical impact on r\,. 

Values of r\, range from sub-parsec to kpc scales. For the 
values of C and / a b s adopted in cases II and III, dark matter 
annihilation could only solve the cooling flow problem for a 
steep (a > 1 .75) profile, in agreement with previous studies 
dTotanill2004t IColafrancesco et al.ll2006l) . On galactic scales, it 
will reduce the fraction of baryons collapsed into stars and pre- 
vent the growth of excessively massive bulges in small systems. 
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Fig. 2. Radius r\, at which heating balances cooling (top), en- 
ergy injection within r\, (middle) and within f2oo (bottom), as a 
function of virial mass. Solid lines correspond to the profile (|4j. 
Results for expression ©, with a varying from 0.25 to 1.75 in 
steps of 0.25, are plotted as dashed lines. 

Ionization would reduce the cooling rate for gas temperatures 
10 4 < T < 10 5 K and increase it for T < 10 4 K. Figure|3]gives an 
estimate of this e ffect by assuming pure thermal bremsstrahlung 
(lEfstathioulll992h . Only the unrealistic Case I would be able to 
completely quench star formation in small (M200 < 10 9 M ) ha- 
los. In all the other cases, r\, < 0.01r2oo. 

Finally, we have tested the dependence on the shape of the 
gas profile by using equation (0) to model the baryonic compo- 
nent. Results depend very weakly on the value of a gas , as long as 
ffgas ^ 0.3. For higher values, which may be attained if cooling 
has already taken place in the outer parts of the halo, dark matter 
annihilation will not be effective in preventing further cooling 
and star formation. Note, however, that this gas will be quickly 
converted into stars, leaving a less dense medium. In fact, such 
mechanism ha s been propo sed to explain the entropy excess in 
galaxy groups (Bryan 2000l). 

4. Conclusions 

Annihilation of dark matt er particles would produce gamma 
ray (see e.g.lBertondl2006l and referenc es therein), X-ray (e.g. 
Bergstrom et al. 2006) and radio (e.g. IColafrancesco & Meld 
200 1 : lColafrancescol2004l) emission. The present study suggests 
that it may also have a noticeable effect on galaxy formation and 
evolution, providing a constant (and relatively powerful) heat 
source at the center of every dark matter halo. 

The magnitude of the effect depends on the physical proper- 
ties of dark matter and its distribution within the halo. Our an- 



2 The time required to ionize all gas within 1%, T mn ~ 
M m {r h ){E i<x ) IE(r b ), where <£ lon > =s 3 x 10 46 erg M" 1 , is of the or- 
der of 1 Myr. 
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Fig. 3. Radius r\, for an ionized gas. Solid lines depict cases I, 
II, and III (see text). Dashed line corresponds to Case II with 
a = 1 .75 in and dotted lines show a = 1 and a — 1 .75 for 
Case III. Dash-dotted line displays the virial radius, r2oo- 



alytical estimates show that, for reasonable values of the model 
parameters, the amount of energy injected into the gas can be 
larger than the cooling rate in the central regions. Only for ex- 
tremely shallow dark matter density profiles or steep gas density 
profiles could all the energy be radiated efficiently. 

Else, cooling and star formation would be completely 
switched off within the radius r\, where heating balances cooling. 
For most of the cases considered, this radius is between 0.01 and 
1 per cent of the virial radius of the object. For the upper limit, 
Case I, our results indicate that no stars could form within 10 per 
cent of the virial radius, and no star at all could form in haloes 
less massive than ~ 10 9 M Q . 

Evaluating the impact of dark matter annihilation on the star 
formation rate outside r\, would only be possible by implement- 
ing dark matter heating in a self-consistent numerical simulation 
of cosmic structure formation. Equation ([1) provides a simple 
prescription to carry out such an experiment and determine the 
maximum amount of heat compatible with current observations. 
Some dark matter scenarios (e.g. Case I) seem to inject too much 
energy for galaxies to form. More realistic models (e.g. cases II 
and III) might actually explain why current models predict many 
more stars than observed, particularly in the central regions of 
galaxies, thus providing an intriguing alternative to more con- 
ventional astrophysical feedback mechanisms. 
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